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ABSTRACT 

We present the first photometric observations of trans-Neptunian objects (TNOs) taken with the 
Kepler space telescope, obtained in the course of the K2 ecliptic survey. Two faint objects have been 
monitored in specifically designed pixel masks that were centered on the stationary points of the 
objects, when their daily motion was the slowest. In the design of the experiment, only the apparent 
path of these objects were retrieved from the detectors, i.e. the costs in terms of Kepler pixels were 
minimized. Because of the faintness of the targets we employ specific reduction techniques and co¬ 
added images. We measure rotational periods and amplitudes in the unfiltered Kepler hand as follows: 
for (278361) 2007 JJ 43 and 2002 GV 31 we get Prot = 12.097h and Prot = 29.2 h while 0.10 and 0.35 
mag for the total amplitudes, respectively. Future space missions, like TESS and PLATO are not well 
suited to this kind of observations. Therefore, we encourage to include the brightest TNOs around 
their stationary points in each observing campaign to exploit this unique capability of the K2 Mission 
- and therefore to provide unbiased rotational, shape and albedo characteristics of many objects. 

Subject headings: methods: observational — techniques: photometric — astrometry — mi¬ 
nor planets, asteroids: general — Kuiper belt objects: individual (2002 GV 31 , 
(278361) 2007 JJ 43 ) 


1. INTRODUCTION 

Kepler has provided incredible results on extraso¬ 
lar planets and planetary systems, as well as on stel- 
lar astrophysics. After analyzing main-belt asteroids 
(|Szab6 et al.l 120151) . here we continue the exploration of 
our own Solar System by pushing the limits of the space¬ 
craft and observing faint, trans-Neptunian objects. 

The Kepler spacecraft is equipped with a 0.95-m 
Schmidt telescope and 42 imaging CCDs, from which 
38 are currently in operation. Kepler was designed to 
be the most precise photometer i n order to detect th e 
transits of numerous small planets (iBorucki et al.l[2M^ . 
It monitored approximately 170 000 targets in the Lyra- 
Cygnus field, at high Ecliptic latitudes, almost continu¬ 
ously for four years during its primary mission. But after 
the failure of two reaction wheels the telescope perma¬ 
nently lost its ability to maintain that attitude. Soon a 
new mission, called K2, was initiated to sa ve the other¬ 
wise healthy and capable space telescope (|Howell et al.l 
I 2 OI 4 II . Since then, Kepler has been observing in shorter, 
75-day-long campaigns along the Ecliptic to balance the 
radiation pressure from the Sun. 

The ecliptic fields include innumerable Solar System 
objects, although most asteroids and planets cross the 
field-of-view quite quickly. TNOs, in contrast, appear 
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to move slowly around their stationary points, and can 
be observed for extended periods of time. The 4"/pixel 
resolution of the CCDs is sufficient to avoid confusion 
even for faint targets. But the observations posed vari¬ 
ous challenges for Kepler. Even the brighter TNOs are at 
the detection limit of the telescope (V = 20 — 22 mag), 
with an expected precision of a few tenths of a mag¬ 
nitude for a single long-cadence (30-min) observation. 
Moreover, TNOs exhibit quadratically increasing proper 
motions away from the stationary points. Since the pixel 
mask allocation is fixed for an entire campaign, the whole 
length of the orbital arc has to be observed continuously. 
Despite these challenges, the K2 mission is in a unique 
position to gather continuous TNO light curves in white 
light. 

The rotational light curve of small bodies in the Solar 
System is determined by the shape and albedo variega¬ 
tions on the surface of the target, and provides infor¬ 
mation on these important characteristics. In addition, 
the collisional history of these bodies is stored in the 
spin state of the objects, probing different eras of the 
Solar Systems’s debris disk evolution depending on the 
size. The largest objects (~ 500 km or larger in diam¬ 
eter) decoupled early from collisional evolution of the 
debris disks, therefore their rotation should reflect the 
cumulative effect of collisions during the formation pe¬ 
riod. Medium size objects (100 — 500 km) likely avoided 
catastrophic impacts with a rotation affected by smaller 
past events, while the smallest bodies (< 100 km) are 
probably remnants of several collisions during the evolu¬ 
tion of the Solar Syst em’s debris disc (for more details, 
see also lLacerdall2005ll . 

In order to study the rotational properties of TNOs 
we proposed a pilot study for K2 Campaign I to observe 
the object 2002 GV 31 (U Ri 22.5 mag) with a small pixel 
mask for a limited time around the stationary point (pro- 
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Fig. 1.— Trajectory of (278361) 2007 JJ 43 and the allocated 661 1 x 11 and 1 x 13 pixel stripes on the CCD module ^17 of the Kepler 
space telescope. In the insets we show three individual exposures along the arc after the stationary point which unfortunately fell off silicon. 
The small inserts show an area of 1.3' x 1.8'. 


posal ID: GO1064). The field of Campaign 2 included 
multiple moderately bright targets from which ultimately 
only (278361) 2007 JJ 43 (V « 20.25 mag) was selected for 
observation (Proposal IDs GO2066 and GO3053). 

2 . OBSERVATIONS 

For 2002 GV 31 a rectangular, 23 x 22-pixel mask was 
allocated by the Guest Observer Office that included 
the object for 16 days around its stationary point, but 
2.8 days were lost to a mid-campaign data download pe¬ 
riod. In contrast, 2007 JJ 43 was covered with a long arc 
constructed from a mosaic of 661 small (1 x 11 or 1 x 13) 
pixel masks, but the region of the stationary point itself 
fell off silicon. This caused a 20.9 d long gap in the ob¬ 
servations, separating the data into a 18.7 d and a 39.2 d 
long section. 

The elongation of the object decreased from 135 to 60 
degrees during the campaign. The public target pixel 
time series files from the Cl and C2 fields were retrieved 
from the MAST archive® for the respective observations. 

Without stabilization around the third axis, the Ke¬ 
pler space telescope slowly rolls about the optical axis. 
This movement causes the field-of-view to rotate slightly 
and it is corrected by the on-board thrusters at every 
5.88h (/corr = 4.08d“^). The spacecraft stores the pix¬ 
els of pre-selected targets only, similarly to the primary 
mission, but larger target pixel masks are allocated for 
less targets in the K2 fields to accommodate the larger 
pointing jitter and less frequent data download periods. 

® https://archive.stsci.edu/k 2 / 


Hence, the data acquisition principles are exactly similar 
to the stellar targets, only the combination of the small 
pixel masks of 2007 JJ 43 required additional steps during 
analysis. 

The attitude of the spacecraft was adjusted by a few 
pixels shortly after the start of both campaigns. This 
did not affect the data of 2002 GV 31 but in the first 49 
frames of Campaign 2, 2007 JJ 43 was too close to the edge 
of the mask and those were discarded. Apart from these 
adjustments, the largest correction movements were al¬ 
ways smaller than a pixel and the root mean squre of 
image shift offsets were 0.3 pixels for both campaigns. A 
closer inspection of the field-of-view movements revealed 
that the direction of motion reversed at the middle of the 
campaign and in some cases the correction maneuver did 
not occur and the telescope drifted for 11.76 h (instead 
of 5.88 h). 

3. DATA REDUCTION 

While in the case of 2002 GV31, the target pixel file 
contained a single stamp of 23 x 22 pixels, the data 
extraction was much complex in the case of 2007 JJ43. 
Here, the total mask related to this object has been split- 
ted to individual target stamps containing only 1 x 11 
or 1 X 13 pixels. Hence, the first step of the reduc¬ 
tion was the reconstruction of the effective section of 
the target CCD areas and the marking of the outlier 
pixels (in order to exclude them from the further data 
processing). In the case of 2002 GV.31 , we employed the 
PyKE package (|Still &: Barclavl[201^ to create individ¬ 
ual FITS frames for each timestamp while the combined 






























TNO observations with K2 


3 



Fig. 2.— Detection of the motion of the faint TNO, 2002 GV 31 
around its stationary point, on 1-day co-added images. These 
stamps are centered at a = 11^28”^38.8®, 5 = -|-05°34^02" and 
cover the 23 x 22 Kepler pixels, yielding a field-of-view of nearly 
1.5' X 1.5'. The image orientation is related to /Cep/er sensors dur¬ 
ing these observations, i.e. north is roughly to the left while east 
is downwards. 


and marked frame series for 2007 JJ 43 was creat ed using 
the utilities of the FITSH^ package ( Pall 1201^ . Prior 
to the light curve analysis (see Sec. 0]), for the subse¬ 
quent data processing steps we employed various tasks 
of the FITSH package (esp. fiign, fitrans, grtrans, 
grmatch, fiphot, fistar). 

Since the target pixel time series do not contain infor¬ 
mation regarding to the time variation of the astrometric 
solutions of the individual smaller stamps, frame regis¬ 
tration had to be based on purely imaging data. Due to 
the fairly strange shape of the field containing the ap¬ 
parent path of 2007 JJ 43 (see. Fig. [T]), automatic source 
identification and cross-matchi ng algorithms bas ed on 
triangulation matching (see also lPal &: Bakosll2006D were 
not found to be effective®. Therefore, we manually se¬ 
lected seven prominent® stars distributed nearly uni¬ 
formly in the field. Namely, we picked four stars located 
in the northern stripe and three stars in the southern 
stripe. Initial pixel coordinates for these stars have been 
acquired using the imexam tool of IRAF^® on adjacent 
frame pairs belonging to events where the pointing of 
Kepler was changed abruptly. We found that this kind 
of manual analysis of 39 images were sufficient in order to 
bootstrap the astrometry and for the cross-identification 
procedure for all of the 3789 individual images. 

After the cross-identification and the derivation of 
the differential transformations were performed, we used 
these transformations to register the images to the same 
system. Then, we stacked every tenth of the images af¬ 
ter correcting of the background variations in order to 
obtain a master frame. This master image was created 
using a median averaging to exclude the effects of the 
cosmic hits, as well as to safely rule out the contribution 
of other moving objects^^. 

Photometry of 2007 JJ 43 was then performed us¬ 
ing differential images (difference between the indi¬ 
vidual images and the aforementioned master image) 
where the aperture centroid coordinates were based 
on Kepler-centric apparent positions retrieved from the 


http://fitsh.szofi.net/ 

® The triangle parameter space is highly biased in this case while 
triangulation matching algorithms works most efficiently when the 
triangles have a nearly uniform distribution in this space. 

® Bright, isolated stars which are close to the center of the stripes. 

IRAF is distributed by the National Optical Astronomy Ob¬ 
servatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under a cooperative agreement 
with the National Science Foundation. 

See http://szofi.net/apal/astro/k2/tno/ for animations. 


TABLE 1 

Photometric data of 2007 JJ 43 and 2002 GV 31 


Object 

Time (JD) 

Magnitude‘s 

Error 

2007JJ 43 

2456893.044947 

20.288 

0.078 

- 

2456893.208405 

20.686 

0.111 

- 

2456893.228837 

20.722 

0.105 

2002 GV 31 

2456839.65 

23.680 

0.606 

- 

2456839.75 

22.862 

0.402 

- 

2456839.85 

22.957 

0.403 

Note. — 

Table is published 

in its entirety 

in the 


electronic edition of the Astrophysical Journal Letters. 
A portion is shown here for guidance regarding its form 
and content. 

^ Magnitudes shown here are transformed to USNO-Bl.O 
R system, see text for further details. 


NASA/Horizons service^^. The instrumental aperture 
was chosen to be relatively small due to the faintness of 
the object. Despite the usage of differential images, the 
residual structures were significant due to the undersam¬ 
pled property of the Kepler camera. Namely, the spline- 
based interpolation throughout the registration process 
always introduces undershootings in case of stellar pro¬ 
files with small FWHMs yielding characteris tic feature s 
on the residual images (see e.g. Fig. 4 in iPall 1^0091) . 
Photometric uncertainties have been estimated consid¬ 
ering photon noise and measured background scatter in 
the aperture annuli. For the photon noise estimation, 
we used the gain values reported in the original FITS 
header, namely 113.3 e“/ADU. As we will discuss later 
on (Sec.|4]), this estimation was found to be reliable. 

In the case of 2002 GV 31 , we followed a similar data 
flow for image reduction, with the exception that the 
frame registration was much easier due to the smaller 
size and the fewer number of stars (Fig. [2]). Namely, the 
number of unambiguously identified background stars 
was three and the mean of the centroid shifts (in both 
directions) was used to obtain the offset values for the 
image registration. In addition, we binned the images 
into intervals of O.ld in order to increase the signal in 
the frames. Hence, each frame on which photometry is 
performed were derived from ^ 5 individual downlinked 
measurements. In all other aspects, the photometry was 
performed identically for the two objects. 

For both objects, photometric magnitud es have been 
trans formed into USNO-Bl.O R system (jMonet et al.l 
[Mil . In the case of 2007 JJ 43 , we found that the 
unbiased residual of the the photometric transforma¬ 
tion between USNO and Kepler unfiltered magnitudes 
was 0.05 mag, while for 2002 GV 31 it was found to be 
0.09 mag. Stars used for both objects have magnitudes 
spanning nearly homogeneously between i? Ri 15 and 
R w 19. Photometric data are displayed in Table |T] 

4. LIGHT CURVE ANALYSIS 

The data produced from the baseline photometry con¬ 
sisted of time instance, magnitude and magnitude errors. 
2224 data points belonged to 2007 JJ 43 , and 129 points to 
2002 GV 31 , respectively. We tested that the photometric 
error estimates by the reduction pipeline are quite re¬ 
alistic. The mean determined photometric error of the 
entire dataset is 0.163 magnitudes, while the mean of the 

http://ssd.jpl.nasa.gov/horizons.cgi 
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Fig. 3.— Fourier transform of the K2 light variation of the 
(278361) 2007 JJ 43 normalized to 1 (red). The frequency axis was 
transformed to the time (period) dimension. In order to demon¬ 
strate that thanks to the motion of the TNO, the pointing ad¬ 
justments were averaged out from the light variation, we plot the 
Fourier transform of the x coordinate of the centroid. The varia¬ 
tions from the TNO rotation and the spacecraft pointing adjust¬ 
ment are clearly separated. 


standard deviation of data points in the boxcar is 0.127 
magnitudes. These values are comparable to each other 
within 30% accuracy. For later processing, the weights 
of the photometric data were set to be proportional to 
1 /cr^, where ai is the formal photometric error of each 
data points. 

We also observed that this approach overweights the 
outliers on the bright side, and may introduce a distor¬ 
tion in the shape of the light curve. This is because error 
sources, such as residuals after image subtraction, scat¬ 
ter toward both directions (too faint and too bright), but 
the mean brightness is biased here, because the formal 
errors of bright-looking objects are much smaller, and 
they are given larger weights. Therefore at the bright 
wing of data distribution, the weights were downscaled 
by a function of the shape 


arctan 


8e-W 


0.3 


( 1 ) 


where d = m — mg, the deviation of the measured m 
magnitudes from the mean value. We selected ^ ss 4 
times the peak of the Fourier spectrum at the rotational 
period (0.15 for 2007 JJ 43 and 0.5 for 2002 GV 31 , respec¬ 
tively). Here the peak of the Gaussian curve is damped 
by the arctan(-) function acting on it, leading to a flat 
plateau around the center of the shape. We are con¬ 
vinced that this rescaling leads to more realistic weights 
than the crude estimates making use of the measured flux 
values. However, in order to avoid biasing the results 
in any way, we compared the results with and without 
weighting, and we experienced that they were practically 
identical - so the modified weights supported merely the 
better visualization. 

The rotation periods were found by a Fourier analysis. 
The Fourier transform of the light curves we re calculated 
by the Period04 code (jLenz fc Bregei][2005l) . To exclude 
the light variation introduced by periodic systematics, we 
also calculated the periodograms of the pointing position 
of the telescope, the actual and the measured subpixel 
positions of the TNOs in each frames. We did not find pe¬ 
riodic signals resulting from the sampling effects (there¬ 


fore subpixel positions), but we identified the peak due 
to the pointing corrections with 4.05 cycle/d frequency 
(Fig.ISj. This is close to the frequency peak belonging to 
a half rotation of 2007 JJ 43 , but they are still well sepa¬ 
rated, proving that the periodograms of the TNOs reflect 
the rotation. 


5. RESULTS 
5.1. (278361) 2007 JJi3 

(278361) 2007 JJ 43 is an outer Kuiper belt object, 
likely in 2 : 1 resonance with Neptune (jMPG 602351 
|20 o 3). Based on its absolute magnitude at discovery 
(iJv = 3.9) and the cluster average geometric albedo 
of p v = 0.13l° °y for objects in the outer Kuiper 
belt (iLacerda et al.l [20l^ . its likely effective diameter 
is Oeff = 610^1™ km. 2007 JJ 43 was also detected by 
the Southern Sky and Galactic Plane Surve y for Bright 
Kuiper Belt Objects (jShepDard et al.l[201 111 , and an R- 
band absolute magnitude of Hr = 3.2 was associated to 
t he object. _ 

iBenecchi fc Sheopardl (|2013f) detected a light curve of 
2007 JJ 43 , with a peak-to-peak amplitude of 0.13±0.02 in 
the Sloan r’band. The most likely light curve (rotation) 
period they found was 6.04 h, or similarly likely its dou¬ 
ble period, 12.09 h. However, other periods (4.83/9.66 h) 
were also possible. The light curve was clearly non- 
sinusoidal. The absolute magnitude in this band was 
found to be Hr’ = 4.17 ± 0.20. 

The period analysis confirme d the most likely solution 
of iBenecchi fc Sheopardl (12013( 1. namely 6.048 ± 0.018 h. 
This is clearly the half of the rotation period, con- 
hrmed by the asymmetry of the phased light curve with 
P = 12.097 ± 0.036 h period (left panel of Fig. S]). 
Weighted means of the magnitude and corresponding 
phases were calculated in a boxcar of 68 data points. 
The standard deviation and the number of points in 
each bins resulted in an expected standard deviation (er¬ 
ror bars) of Ri 0.015 magnitudes for the binned points. 
The smoothed light curve shows two humps with equiv¬ 
alent minima and slightly different peaks. The full am¬ 
plitude is 0.100 ± 0.005 magnitude, the difference be¬ 
tween the minima is 0.02 magnitudes. The asymmetry 
and the non-sinusoidal shape of the light curve reflects 
the uneven surface structures. These features are quite 
common for Solar System objects, e.g. a very similar 
light curve shape is observed for the main belt aster- 
oid (52) Europa at certa in observing circumstances (e.g. 
iMichalowski et aT1l2004[l . We can also conclude that the 
presented K2 measurement and the previous photometry 
of 2007 JJ 43 in IBenecchi fc Sheopardl (|2013H shows sim¬ 
ilar characteristics since in the case of TNOs, the phase 
and aspect angles varies rather slowly on the timescale 
of these observations. 

5.2. 2002 GVgi 

2002 GV 31 is a dynamically cold classical Kuiper belt 
object. A recent evaluation of Johnson R-band MFC 
data and a conversion using V — R = 0.59 ± 0.15 (aver¬ 
age for classicals) resulted in a V-band absolute mag¬ 
nitude of Hv = 6.1 ± 0.6. 2002 GV 31 was observed, 
but not detected in any bands of the PAGS camera on 
board the Herschel Space Observatory, in the frame¬ 
work of the ’’TNOs are Cool!” Open Time Key Pro- 
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Fig. 4.— Left: phased light curve of (278361) 2007 JJ 43 (red points), 
light cu rve of 2002 GV.s i • _ 

gram ijMiiller et al.l]200^ . Using the upper limits set 
by the non-detections, a modelling of the thermal emis¬ 
sion of the target resulted in an upper limit of its size 
{D < 180 km), an d a lower limit on its geom etric albedo 
(0.19 < pv, see iVilenius et al.l I 2012 L ^0141) . No light 
curve information is available for this target. 

The first light curve of this object is shown in the right 
panel of Fig. IH phased with the most likely rotation pe¬ 
riod of 29.2 ± 1.1 h and showing a double peak. Since we 
had 129 photometric points from K2 in this case, they 
were binned into 13 points during one rotation, and the 
error bar of the binned points resulted to be 0.06 mag¬ 
nitudes. The full amplitude of the light variation is 0.35 
magnitudes, which can be a sign of a more exposed as- 
phericity than for 2007 JJ43. The most interesting fea¬ 
ture of this TNO is its long rotation period, as most 
TNOs and Centaurs have rotation periods shorter the n 
a day (|Duffard et al.l[200^ iBenecchi fc SheDDardll201^ . 

Longer periods are expected from synchronously 
locked, tidally e volved binary systern s (see Pluto/Charon 
or Sila/Nunam, iGrundv et al.ll2012D. It can explai n the 
slow rotation of 2010 WGq (iRabinowitz et al.ll20ll) and 
can be plausible for 2002 GV 31 either. Moreover, ob¬ 
servational statistics can explain binarity more likely for 
dynamically cold TNOs (|Noll et al.ll20f)8[ ). Ii i addition, 
statis tics seem to prefer close orbits (Fig. 4 in iNoll et alJ 
I 2 OO 8 D which can also yield mutual events and hence 
strong light curve variations. 

6 . CONCLUSIONS 

We clearly detected the rotational signal of two TNOs 
with Kepler during the K2 extended ecliptic survey mis¬ 
sion. Since these objects are moving through the field- 
of-view, photometric data reduction requires special care 
in order to retrieve accurate flux information. Although 
the image processing needs complex steps, our results 
clearly show that the moving nature of these targets 
are also advantageous. Namely, differential photometric 
techniques strenuously reduce the effects of background 
sources, even in the crowdness of Kepler fields and in 
the case of lower imaging resolution and undersampled 
stellar profiles. In addition, the employment of moving 
apertures removes the frequency aliases caused by the 


The bold symbols with the error bars are binned values. Right: K2 

pointing corrections of the spacecraft in every ^ 6 hours. 

All in all, we successfully demonstrated that it is feasi¬ 
ble to observe faint TNOs around their stationary points 
with Kepler/K2 and this mission provides a unique way 
to monitor them continuously for many weeks. The sci¬ 
entific benefit of observing TNOs is the opportunity to 
provide an unbiased sample of rotational light curves - 
and hence data about shape, albedo and surface charac¬ 
teristics. This kind of information aids us to understand 
the nature and evolution of the outer Solar System. 

Future space missions, like TESS and PLATO are not 
well suited to this kind of observations, therefore we en¬ 
courage to include the brightest TNOs in each observing 
campaign to exploit this unique capability of the K2 Mis¬ 
sion. 
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